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Abstract 



Considering the constraints from the experimental data on // — > ey, — > 
3e, /i — e conversion etc., we analyze the Lepton Flavor Violating decays 
4>( J/ty, T(lS')) — > e + /i _ (yU + r _ ) in the scenarios of the minimal supersym- 
Q-i metric extensions of Standard Model with seesaw Mechanism. Numerically, 

r-| I there is parameter space that the LFV processes of J/\I/(T) — > h + t~ can 

reach the upper experimental bounds, meanwhile the theoretical predictions 
on n — > ej, n — V 3e, \i — e conversion satisfy the present experimental 
bounds. For searching of new physics, Lepton Flavor Violating processes 
J/^Y) — > fi + r~ may be more promising and effective channels. 
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1. Introduction 

As an evidence to discover new physics beyond the Standard Model (SM), 
searching for Lepton Flavor Violating (LFV) processes in charged lepton 
sector have attracted a great deal of attention. The theoretical predictions 
on those lepton flavor violating processes are suppressed by small masses 
of neutrinos in SM, and exceed the detecting extent of experiment in near 
future. Nevertheless, the corrections to the branching ratios of LFV decays 
cf) — > e + fi~, J — > fi + T~ and T — > fi + r~ are enhanced by the new sources 
of LFV in various extensions of the SM, such as grand unified models 
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supersymmetric models with and without R-parity 0, left-right symmetry 
models 0] etc. Although nonzero neutrino masses supported by the neutrino 
oscillation experiments j4| imply the non-conservation of lepton flavor, it is 
very important to directly search the LFV processes of charged lepton sector 
in colliders running now. 

Using a sample of 5.8 x 10 7 J '/\& events collected with the BESII detector, 
Ref.0] obtains the upper limits on BR( J/W y.r) < 2.0 x 10" 6 and BR(T 
yur) < 8.3 x at the 90% confidence level (C.L.). Adopting the data 
collected with the CLEO III detector, the authors of Ref. |6| estimate the 
upper limits on BR (T(15) -> fir) < 6.0 x 10" 6 , BR (T(2S) ->■ \ir) < 
1.4 x 10" 5 and BR (T(3£) ->■ fir) < 2.0 x 10~ 5 respectively at the 95% C.L. 
Additionally, the study of LFV processes involving light unflavored meson is 
an effective way maybe to search for new physics beyond the SM, and the 
SND Collaboration at the BINP (Novosibirk) presents an upper limit on the 
4> — > e + ji~ branching fraction of BR (0 — > e + fi~) < 2 x 10~ 6 0. 

In literature, several stringent limits on LFV decays of both light and 
heavy unflavored mesons are derived already. Assuming that a vector boson 
Mi (Mi could be either a fundamental state, like the Z , or a quark-antiquark 
bound state like the 0, J/^f, T) couples to /i T e ± and e T e ± as: 

C eff = g MiPe Wn eM ? + 9 Mi e^ eM i + h - c - > (!) 

where g M ee and g M , e denote the corresponding couplings of a meson to lepton 
flavor conservation and violation currents, and by unitarity its exchange con- 
tributes to fj, — >■ 3e, the authors of Ref. jsj deduce upper bounds on the LFV 
decay of mesons from the LFV process fi — > 3e. Under a similar assumption 
that a vector meson Mj couples to yU =F e ± and NN as: 

£ e// = (Cvn^ + tTn^Wt + 9mnn n^nm? + h.c. , (2) 

where N is a nucleon, ^y A are effective vector and axial couplings of a meson 
to the LFV lepton currents, authors of Ref. H studies the LFV decays of 
vector mesons by taking account of the experimental constraint on \x — e 
conversion. It shows the constraint from \x — e conversion on LFV decays 
of vector mesons is more stronger. Likewise, authors of Ref. jsj also deduce 
upper bounds on other LFV decay of mesons from the LFV processes r — > 3e 
and r — y 3/z. Making the assumption that fermion mixing and mass hierarchy 



originate from mass matrix rotation, authors of Ref. [lOj also get some upper 



limits on the LFV decays of heavy unflavored mesons and Z boson. Searching 
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for new physics beyond the SM is also a goal of LHC In LHC, vector mesons 



can be produced by photo fusion 11 



In SM,the LFV decays mainly originate from the charged current with 
the mixing among three lepton generations. The fields of the flavor neutrinos 
in charged current weak interaction Lagrangian are combinations of three 
massive neutrinos: 

£ = V T2(x)w lL (x)W»(x) + h.c, 

^ 2 i= 



^2 (UmnsJ (3) 



where g 2 denotes the coupling constant of gauge group SU(2), uil are fields 
of the flavor neutrinos, v iL are fields of massive neutrinos, and Umns corre- 



sponds to the MNS neutrino mixing matrix [L2j,ll3|].In the standard parametriza- 



tion |l4|], the leptonic mixing matrix is given by: 

cic 3 C3S1 
U M ns = I -c\s z s 2 e %s - c 2 si c x c 2 - s 1 s 2 s 3 e lS 
S1S2 - CiS 3 c 2 e J<5 cis 2 - sic 2 sie t5 




xdiagU^/ 2 ,!,^^ 2 ), (4) 



where s(c)i = sin(cos)6 ) i 2 , s(c) 2 = sin(cos)6 , 23 , s(c) 3 = sinfcos)^. The phase 
5 is the Dirac CP phase, and $j are the Majorana phases. A global fit of the 
neutrino oscillation data points out: Q\ 2 ~ 34° and # 2 3 ~ 45°. Recently, the 
observing v e disappearance in reactor experiments Daya Bay [l5[ and RENO 



16( have definitely established that #13 > at ~ 5a level. The Daya Bay 
and RENO have measured sin# 13 ~ 0.024 and sin^ 13 ~ 0.029 , respectively. 
However, the theoretical predictions on branching ratios of any LFV decays 
are suppressed strongly by the tiny neutrino masses in SM and fall out the 
reach of experiment in near future. In this work, we analyze the LFV decays 
: — > e + /i~, J/'fy — > fi + T~ and T(IS') — > fi + r~ in the framework of ac- 
comodating supersymmetry with type I seesaw mechanism simultaneously. 
With the accumulation of events on BEPC [5( and SuperKEKB |17| . the 
updated experimental data on those LFV decays maybe constrain the con- 
cerned models more stringent. To shorten the length of text, we just present 
the upper bounds on those branching ratios of p(u, J/^, T) — > e + /i~ under 
our assumptions on parameter space. 
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The paper is organized as follows. In Sectionf2J we firstly provide a simple 
overview for the origin of lepton flavor changing and corresponding interac- 
tion lagrangian in the framework of MSSM with type I seesaw mechanism. 
Then, as an example,we derive the analytic results of amplitude for one di- 
agram in detail. The numerical results are presented in Sectionj3l and the 
conclusion is drawn in SectionJH All the simplified amplitudes corresponding 



to the Feynman diagrams in Figfj] and FigfS] are given in |Appendix A 



2. Formalism 

In the minimal supersymmetric extension of SM with R-parity conserva- 
tion, the general form of the superpotential involving the lepton and Higgs 



superfields is given by [18 



Wmssm = e i '{^HlH 2 j +Y l IJ H 1 i L I j R J ), (5) 

where /i is the mu-parameter, the 3x3 matrix Y\ is the charged lepton 
Yukawa couplings. For convenience, we assume Y L IJ = Y/S 1,7 (I,J=1,2,3) in 
this work. Then, the relevant soft supersymmetry breaking terms involving 
the slepton sector and sneutrino sector are: 

-A , f J H?*L I i R J - h.c. , (6) 

where m 2 L is left 3x3 soft slepton mass matrix, m 2 R is right 3x3 soft slepton 
mass matrix, the 3x3 matrix Ai is the trilinear scalar couplings, the 3x3 
matrix A[ is the non-standard trilinear scalar couplings, respectively. The 
LFV interactions mainly originate from the potential misalignment between 
the leptons and sleptons mass matrices in the MSSM. In other words, the 
sources of LFV are the off-diagonal entries of the 3x3 soft supersymmetry 
breaking matrices mf,, m 2 R , Ai and A[ in 6 x 6 slepton mass matrix Mj?, which 
are listed below: 

m a" = w » + mi,., + {mlY , , (7) 

J ll 8sici 2 

Ml)" = -?M^l s u + mi s » + (m y» , ( 8 ) 



RR Acl 

(m£) U lr = ±{^^6" - A!{ J ) + v x Al J ) , (9) 
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where s w =sin 6 l w ,c w =cos 6 W with 6 W denoting the Weinberg angle, and t?i j2 
are the non zero vacuum expectation values (VEVs) of two Higgs doublets. 
In the minimal supersymmetric extension of the seesaw extended SM 



191 . |2Q| . |21| . |22L 1231 . [24|, there are three generation right handed neutrino 
superfields N 1 (I = 1, 2, 3) with zero hypercharge. The most general 
form of the superpotential involving the lepton and Higgs superfields in the 
R-parity conserving scenario is given by: 

W = W MSS M + e ij Yj J H?LiN J + Im ij N i N j , (10) 

where Y v is the 3x3 neutral lepton Yukawa coupling, M is the 3x3 Majarana 



mass matrix. Here, we adopt the parameterization in [19| to reproduce the 
PMNS mixing matrix: 

3 

{Y v f = 1 V2(m ul M u ^ 2 R jk (Ut INS ) tk /v 2 , (11) 

k=l 



where Umns is the MNS mixing matrix in Eq.(j4]), m^(z = e,/i, r) are the 
masses of left handed neutrinos, and M^(i = e, /i, r) are right handed neu- 
trino masses. Furthermore, R is an arbitrary orthogonal matrix [25[ deter- 
mined by three angles a±, a 2 , a 3 : 

(C2C3 -C1S3 - S1S2C3 S1S3 - C1S2C3 \ 
C2S3 cic 3 - S1S2S3 -S1C3 - cis 2 s 3 , (12) 
S2 SiC 2 CiC 2 / 

in which q = cos and Sj = sin aj, i = 1, 2, 3. In the scenarios of MSSM with 
Seesaw mechanism, the corrections from right handed Majorana neutrinos to 
the branching ratios of vector meson LFV decays can be ignored since they 
are suppressed by the huge masses of right handed neutrinos. In addition, 
the mass term for the light sneutrinos is given by: 

-£T ss = \{~vl vL*)Ml(*iy (13) 
with /, J = 1, 2, 3 are the indices of generation, and the 6x6 mass matrix 



is 



. r I (M 2 LC y j (M 2 L * v y j , 

M ° ~ 1 {Ml v ) IJ (M? c )" J ■ 
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Here M\ c and M| y are 3x3 matrices. If M\ v = 0, the six light sneutrinos 
are comprised of three sneutrino-antisneutrino pairs. If M\ v ^ 0, the lepton 
number is violated and the sneutrinos and antisneutrinos can mix and yield 
six non-degenerate sneutrinos. The elements of M\ c and M\ v are given by, 
in a simple form at GUT scale: 

(M 2 LC ) IJ = (m 2 L ) IJ + l -Ml cos (15) 

(Mi v y j = _ MVcot/3 (nM . lyJ)JJ; g) 



where M is the right handed neutrino mass matrix in Eq. ffl0|) and tan (3 = 
V2I V\. M\ c reproduces the well known 3x3 light sneutrino matrix in MSSM. 
In the CP-base 

*F^(M), (17) 

4-" = -^(s[-f£-), (18) 

the mass term for the light sneutrinos is rewritten as: 

1 



r,(+) J 



= ± ( ^)Ml[ ^ ) . (19) 



Here, the 6x6 mass-squared matrix Ail is 



The effective squared-mass matrix can be diagonalized by 6 x 6 orthogonal 
matrix, Z# via: 

^=K,<... I <) I (21) 

where Si(i = 1, ...,6) correspond to the physical sneutrino mass eigenstates. 
The sneutrino interaction eigenstates, v 1 , can be expressed in terms of the 
physical sneutrino mass eigenstates Sk by: 

^ = 4£(^ fc + ^V*, (22) 

V 2 fc=i 
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Figure 1: The penguin diagrams of the LFV processes cj>(J/^f,T) — > e + fi ) in 

MSSM with seesaw mechanism. 



Correspondingly, the relevant Lagrangian is given as: 



C 



+ 



{~7^ (Z]/s w + Zn°w) + z[ J+3) Z%fjP L 



V2s 
— e\/2 (j + 3)i 



-Xf (f + tfi^P*) (z? - iZ^yS, + /i.e., (23) 

where Z± are the mixing matrices of chargino sector, Zjv, Z^ and Zy are 
the mixing matrices of neutralino sector, slepton sector and sneutrino sector, 
respectively, e 1 denote the SM charged leptons. Lf and Sj denote the 
sleptons and sneutrinos. x? an d X*j stand for the charginos and neutralinos. 
Pl/r — o(l -¥ 75)- The relevant Feynman diagrams contributing to the LFV 
decays are presented in Figfj] and FigfJJ 

In the quark picture, mesons are composed of a quark and an anti-quark. 
As we analyze those LFV processes mentioned above, we do not want to 
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(a) (6) 

Figure 2: The box diagrams of the LFV processes ^(J/'S, Y) — > e + /i~(/i + r~) in the MSSM 
with seesaw mechanism. In concrete calculation, the corrections from cross diagrams 
should be included also. 



calculate the complicated loop integrations at quark-gluon level since the lack 
of a completely reliable way to calculate the non-perturbative QCD effects. 
We adopt a phenomenological model where the amplitude of hard process 
involving a s-wave meson can be described by the matrix elements of gauge 
invariant nonlocal operators, which are sandwiched between the vacuum and 



the meson states. For our case, the matrix is given by [26 



(0\q(y)r[y,x}q(x)\(f>) , (24) 

here T = 7 M or 7^75 is a generic Dirac matrix structure, x and y are the 
coordinates of quark and anti-quark. The distribution amplitude of vector 
meson in leading-order is defined through the correlation function : 







muffin m 



+^' v ' fUz^'Pv' -e*wyW(u) > (25) 

Z \ /J pa 



J /3a 

where the momentum of is on-shell, i.e. p 2 = m^, is the polarization 
vector, and /J are the decay constants of <fi meson, 0y and are the 
leading-twist distribution functions corresponding to the longitudinally and 
transversely polarized meson, respectively. For the cases of J/\& and T, there 
are similar distribution amplitudes. The integration variable u corresponds 
to the momentum fraction carried by the quark, u = 1 — u stands for the 
momentum fraction carried by the anti-quark, a and are the indices of 
matrix elements, and Nc is the number of colors, separatively. Since the 
leading-twist light-cone distribution amplitudes of meson are close to their 



asymptotic form 27], so we set 0ii = 4>± = Mu) = 6u(l — u) 
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Figure 3: Take Figflja) for example, and particle momentum is fixed. 



Taking the diagram in FigJT] as an example, we show how to write the 
corresponding correction to the LFV decay <fi — > e + /i~ in MSSM with seesaw 
mechanism. At quark level, the relevant amplitude is written as: 



—e' 



9iiv 



d D k {zfzf 



444 J (2ix) D [fa+nY-mlW-m*] 



x 



v s {p2)Y{Pl - |4K(gi) 

[(ps - k) 2 - m)][(p 4 + k) 2 - m 2 } 



Z L l Z iV S W + C W) + *l L N 



+ (ztZ 1 N k + Y 2 Z 2 ^Z^)p L ) ( J + m x o) 



P 



+ {z^Z^+YlZ^Zf*)p R }v e {p,). 



Pl 



(26) 



In the frame of center of mass, one can write down the amplitude at hadron 
level using Eq. fl25l) : 



A 



H 



d D k 



' Z 2i Z^* 



24^) 



24iV c 44 J (2n) D [(pi + p 2 ) 2 - m|] [A; 2 - m 2 ] 
Um^Asl - 3) e(p) • (p 3 + p 4 ) _ 

[(p 3 - jfe)2 _ m 2][(p 4 + fc)2 _ m 2] 



x { \zf*{Z lk *s w + Z 2 N k *c w ) + Y 2 ZtZfj 
+ (ztZ 1 N k + Y 2 Z 2 -Zf)p L \ ( ]k + m x o 



P 



X • 



Z 1 l{Z 1 n s w + Z'nCw) + Y^Z^ Z 3k 



N 



Pr 



+ [Z^ j Z 1 N k * + Y, l Z) j Z ik *)P R \v e (p 3 ). 



(27) 



Applying the high energy physics package FeynCalc [28j, one can simplify 
the amplitude in terms of invariant Passarino-Veltman integrals [29| : 



A 



H 



ie 2 Ti 2 f <t> m^Asl - 3) 
24JV r «2 cl (ml - ml 



6,6,4 



i,j,k=l 



X 



AguM{c ime {A*A* *P L + A%Af*P R ) 



+C 2 



im, 



A k s j Af* - m„A k jAf*)P L 



+ (m e A k Mf*-m fl A^Af*)P B 



m 



+ M^ ! r - lAl 3 AT)P L MPS) 



(28) 



All of integrals can be calculated through another high energy physics pack- 
age LoopTools 30] . In a similar way, we can write down the corrections from 
other diagrams in FigJH and Figj2] at hadron level, and list the simplified 
amplitudes in Appendix A. 

Using the summation formula 



A=±1,0 



ni 



2 ' 



(29) 



we express the branching ratio of <fi — > e + /i as 

^[mj - (m e + m M ) 2 ][mJ - (m e - m^) 2 ] 



Br((J) — > e 



16™^ 



x 



(30) 



in which is the total decay width, Ai are the amplitudes in Appendix A. 
The branching ratios for J/^/(T) — > fi + r~ can be formulated in a similar 

way 
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3. Numerical Analysis 

In the numerical analysis, we adopt the following value for mass of mesons 
M = 1.019GeV, Mj/* = 3.096GeV, M T = 9.460GeV. For the decay con- 
stants, we take fa = 0.231GeV, f J/9 = 0.405GeV, f T = 0.715GeV 0. 
Furthermore, the electromagnetic coupling is determined by a(m z ) = 1/127. 
Coinciding with the neutrino oscillation data and not losing generality, we 
always assume the lightest neutrino mass as: m u e = 1.0 x 10^ 14 GeV, and 
the masses of three neutrinos satisfy following relations from experiment: 
Am^ = 8.0 x 10~ 5 eV 2 , AmJ tra = 3.0 x 10~ 3 eV 2 . Here, we also assume three 
right handed neutrinos are degenerate, i.e., M„? ~ M„m ~ M v t ~ M , M is 

R R R 

the mass scale of three right handed neutrinos. The recent results of the LHC 
experiments indicate that the lower limit of the squark mass is roughly given 



as 800 GeV 32J. Not losing generality, we assume the degenerate spectrum in 
scalar quark sector (m^u = (mfj)jj = (m 2 D ) I j = rrASjj = ITeV 2 , A 1 ^ = 
(J, J — 1, 2, 3) at GUT scale to satisfy the constraint. Through the calcu- 
lation of mass spectrum and mixin g m atrices, a publicly available fortran77 
program SUSY_FLAVOR is used 3^ • 



In our numerical analysis, we assume that the gaugino masses are GUT- 
related, that is, 

M 1 = ^4m 2 , M 2 = ^M 3 w Im 3 . (31) 
3c 2 a s 3 

In order to decrease the number of free parameters involved in our calculation, 
we suppose that the diagonal entries of two 3x3 matrices m 2 L , m 2 R in Eq.flH]) 
are equal {m 2 L )n = (m 2 R )jj = m|, where 1 = 1, 2, 3. Now, the only sources 
of LFV are off-diagonal entries of the soft breaking terms m 2 , m 2 R and Ai. 
Those off-diagonal entries of 3 x 3 matrices to 2 , m 2 R are parameterized by 



mass insertions as in 33 



i.i 



<) = *WK)W J . ( 32 ) 

IJ 



m R 



^y(m 2 R ) n (miy j . (33) 
isume 

soft breaking coupling is parameterized by 



where /, J = 1, 2, 3. We also assume 5f J = 5 R . Meanwhile, the trilinear 



A 1 / = ajY/ { (m 2 L Y^m 2 R y^ (34) 



A[ J = S&j2{ml)"(rr%)» (35) 
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At first, we discuss the LFV decays of vector mesons <p — > e + fi~ . The correc- 
tions from Higgs to the LFV branching ratios of vector mesons <fi — > e + /i~ can 
be neglected safely since they are suppressed by the tiny masses of leptons. 

In the MSSM with type I seesaw, the LFV processes originate from the 
mass insertions 8 l [, 5%. The most challenging experimental prospects arise 
for the CK(fi — e) in heavy nuclei such as titanium (fjTi). The experimental 
upper bounds on the conversion rate reach CR (/i — e, Ti) < 4.3 x 10 



-12 



14j. In the MSSM with type I seesaw, the conversion rate in nuclei can be 



calculated by [20 



T(fi + X -» e + X) 
r(/i + X — > capture) 



4a 5 ^|F(g)|W Z (A^ - A%) 



+ 



-{2Z + N)D^ -{Z + 2N)Dl 
Z(A? - A%) - (2Z + N)D* 



-{Z + 2N)D 



R 



2-i 



(36) 



where Z and N denote the proton and neutron numbers in a nucleus, F(q 2 ) 
is the nuclear form factor and Z e ff is an effective atomic charge, A\^ stand 



for the contributions from penguin-type diagram, D^f stand for the contri 



u.d 



butions from box- type diagrams. In ^Ti, F(q 2 ) ~ 0.54 and Z, 



eff 



34j . After a scan over the parameter space, we will discuss 



e conversion for two cases: (I) 8 12 dominance, 8 2 ^ 



5? 



0; (II) 5?5_ 



17.6 
and 

23^13 



dominance, 8 12 
aj — 1 and 5jf R 



■ 0. We assume tanf3 
0. 

ITeV, M n 



10, fj. = 200GeV, M 2 = 200GeV, 

Case (I) Taking m & = ITeV, M = 10 10 GeV, we plot the theoretical 
prediction of BR (0 — > e + /i~) (solid line) and CR^+^Ti — > e+^Ti) (dash 
line) versus 8 12 in FigJU^a), where the gray shadow is the excluded region for 
CR(/i +i Ti ^ e +f 2 Ti) . The CR(/i +| Ti ^ e +f 2 Ti) exceeds the current 
experiment limit at 8 12 ~ 2.0 x 10~ 6 . The parameter space of 8 12 has been 
highly suppressed with respect to the prediction of CR(/i + 2 % Ti — > e + 2 % Ti). 
Both BR (<j) e + /i~) and CR(/i +f 2 Ti ^ e +f§ Ti) tend to be not sensitive 
to S] 2 when its value is below 10~ 7 . In 20(, the authors investigate the LFV 
processes /i —> ej, fi — > 3e and deduce a constraint with 5 12 < 3 x 10~ 4 . 
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(c) MJGeV] 

Figure 4: Case (I):BR (<j> -> e + fj,~) (solid line) and CR(/z +|| Ti -> e +f| Ti) (dash line) 
vs mass insertion ^j 2 , slepton mass sector m ^ and right handed neutrinos mass scale Mo . 
The shadow is the excluded region for CR(/x +£| Ti ^ e Ti). 
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(c) MJGeV] 

Figure 5: Case (II):BR (0 -> e+/i") (solid line) and CR(^ +i Ti -> e +|| Ti) (dash line) 
vs mass insertion (J^ 3 ^ 3 , slepton mass sector rag. and right handed neutrinos mass scale 
M . The shadow is the excluded region for CR(/x +22 Ti ^ e Ti). 
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Comparing with the constraint on from /x — > ej, \i —> 3e, It displays the 
constraint from /i — e conversion is more stronger. 

Taking M = 10 10 GeV, = 10~ 6 , we plot the theoretical prediction 
of BR (0 -»■ e + fi~) (solid line) and CR(> +f| Ti -»■ e +f| Ti) (dash line) 
versus slepton mass sector in Fig|4jb), where the gray shadow is the 
excluded region for CR(/i +22 Ti — > e +22 Ti) . Different to Figfjja) and 
Figgjc), Both BR (0 -> e + /i~) and CR(/i +f 2 Ti ^ e +*J Ti) decrease as 
varies from 0.5TeV to 3TeV. For lower slepton mass, the prediction on 
CR(/i +22 Ti — > e +22 Ti) is also out of the current experiment limit. 

Taking = ITeV, 8)? = 10~ 6 , we plot the theoretical prediction of BR 
(0 — > e + (solid line) and CR(/i +22 Ti — > e +22 Ti) (dash line) versus the 
right handed neutrino mass scale M in FigfJJc), where the gray shadow is 
the excluded region for CR(/i +fj Ti -> e +f 2 Ti). Both BR (0 e + n~) 
and CR(/i Ti — > e +22 Ti) also show a strong dependence on M in range 
of Mo > 10 10 GeV, but most part is ruled out by the constraint from [i — e 
conversion. When M < 10 10 GeV, the dependence becomes weaker and 
weaker. 

Case (II) Taking = ITeV, M = 10 10 GeV, we plot the theoretical 
prediction of BR (0 -> e+//~) (solid line) and CR(/i +i Ti ^ e +| Ti) 
(dash line) versus 5| 3 5i 3 in FigJSfa), where the gray shadow is the excluded 
region for CR(/i +22 Ti — > e +22 Ti) . There is a sharp decrease around 
^zMi 3 ~ 117 x 1CT 6 with a minimum BR (0 — > e + /i~) of order about 10~ 23 , 
which is about two orders smaller than the most stringent prediction in 



In 20|, the authors also give a expected value for 8j?5}? deduced from the 
processes \i — > and fi — >■ 3e, which is (~ 1CT 6 ) and compatible with ours. 

Taking M = 10 10 GeV, 8f8f = 1.2 x 1(T 6 , we plot the theoretical pre- 
diction of BR (0 e + /i") (solid line) and CR(/x +| Ti ^ e +f| Ti) (dash 
line) versus in FigJS^b), where the gray shadow is the excluded region for 
CR(/i +i Ti -» e +i Ti). Here and followed, we will assume 8f = 4 x 10" 3 , 
(5f = 3 x 10 -4 and these are also expected values for 8 2 ^ and 8^ reported 
in [io[ evaluated from LFV decays r — > ^7 and r — )■ e7- We also find a res- 
onating absorption around vn^ = ITeV that originates from the interference 
between the corrections from sneutrino sector and that from charged slepton 
sector. Fig. 5(b) displays that no constraint on mg has arisen with respect 
to fi — e conversion and LFV decay — > e + yT . 

Taking = ITeV, 8f^8}^ = 1.2 x 10~ 6 , we plot the theoretical prediction 
of BR (0 e+/i-) (solid line) and GR(ji +f 2 Ti ^ e +f 2 Ti) (dash line) 
versus M in Figfj^c), where the gray shadow is the excluded region for 
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CR(/i +22 Ti — > e +22 Ti). There is a resonating absorption around M = 
10 10 GeV that originates from the interference between the corrections from 
sneutrino sector. 

Comparing Case (I) with Case (II), we find: (i]In Case (II), the prediction 
for BR (0 e + /i _ ) are more compatible with Q. In jsj and jsj, the con- 
straints are BR (aj> -> e + fx~) < 4.0 x 10" 17 and BR (0 -»■ e+//~) < 13 x 10_21 - 
It is noted worthwhile that the prediction in |8j also satisfies the constraint 
from fx — e conversion, even if it is derived by the constraint from fx — > 3e. 
(ii) Compared with 5f (~ 10~ 3 ) and 5f (~ HT 4 ), the value for 8f (10~ 6 or 
little) is so small that it can be neglected. 

Then, we will investigate meson decays J/\I/(T) — > li + t~ in Case (II) not 
only for reasons above, but also for the aim to generate a large enough BR 
(J/\I/(T) — > fx + r~) to be observed in experiment. As it displays in Figj6] 
and Figd the mass insertion 5£ 3 affects the theoretical evaluation of BR 



(J/\I/(Y) —7- h + t ) strongly. In formula, there is a simple relation 20 



BR(r^3fx) a /16, m T 14 \ n nnn , OPT . 

) 4- ~ — [—In — ~ 0.003. (37) 

BR{t Lii) 8vr V 3 2m M 9 / y ' 

So, we just consider the constraint from r — > fxj. 

(I) J/* -> /i+r" 

Taking mg = ITeV, M = 10 10 GeV, we plot the theoretical prediction 
of BR (J/ty —> /U + r~) (solid line) and BR(r — > fxj) (dash line) versus 5]^ in 
FigEfa). Taking M = 10 10 GeV, 8f = 4 x 10" 3 , 5i 3 = 3 x 10" 4 , we plot the 
theoretical prediction of BR (J/\I/ — > /U + r~) (solid line) and BR(r — > fij) 
(dash line) versus in Figj6]^b). The gray shadow is the excluded region 
for BR(r — > fxj) from the experiment. A linear relationship is displayed 
between BR (J/^f — > /U + r~) and <5f 3 in logarithmic scale, which shows the 
great dependence of BR (J/^f — >■ li + t~) on 5| 3 . We also investigate the 
dependence of BR (J/^ — > /i + T~) on M and it shows BR (J/^ — > /i + r~) is 
not sensitive to M . 

In Q and [13], the authors calculate the BR («//\& — > /^ + t~) with BR 
(J/* ->■ /U+r-) < 4.1 x 10" 9 and BR (J/tf -> /i+r") < 6.3 x lO" 10 , that 
is three or four orders below the experiment limit 2.0 x 10~ 6 . Under our 
assumption, m^ = ITeV, M = 10 10 GeV, 8f = 4 x 10" 3 , <5i 3 = 3 x 10" 4 , we 
get the BR (J/^ — > /U + r~) can be enhanced as large as 1.6 x 10~ 7 , that is 
more promising to detect directly in experiment in near future. 

(II) T -)• ix+t- 
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1.5 2.0 
(b) MJTeV] 



Figure 6: BR (J/* — > h + t ) (solid line) and BR(t — > ^7) (dash line) vs mass insertion 
^f 3 , slepton sector m^. The gray shadow is the excluded region for BR(r — > ^7). 
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(b) M e [TeV] 

Figure 7: BR(T — > /i + T~) (solid line) and BR(r — > ^7) (dash line) vs mass insertion Jf 3 , 
slepton sector mg. The gray shadow is the excluded region for BR(r — > /ry). 
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Table 1: The upper bounds on the branching ratios of vector bosons 



Decay 


Experiment 


Ref.[8] 




Ref.[9] 




Our prediction 


<f) — > e + fi~ 


< 2.0 x KT 6 


< 4.0 x 10' 


-17 


< 1.3 x 10- 


-21 


< 5.0 x 10~ 20 


p —7- e + fi~ 




< 3.8 x 10" 


-20 


< 3.5 x 10- 


-24 


< 1.0 x 10~ 20 


uj — y e + /i~ 




< 8.1 x 10" 


-16 


< 6.2 x 10- 


-27 


< 1.8 x 10~ 20 


J/^ -)■ e + /i~ 


< 1.1 x 10~ 6 


< 4.0 x 10" 


-13 


< 3.5 x 10- 


-13 


< 1.9 x 10~ 18 


T e+/i" 




< 2.0 x 10" 


-9 


< 3.9 x 10- 


-6 


< 3.6 x 10~ 18 


J/\& ->■ /i + r~ 


< 2.0 x 10~ 6 


< 4.1 x 10- 


-9 






< 1.6 x 10~ 7 


T ->• h+t- 


< 6.0 x 10~ 6 


< 7.2 x 10- 


-5 






< 5.3 x 10~ 7 



Similar to the process J/^f — > h + t~ , the LFV decay T — y /i + r~ has the 
same behavior as a function of (Jf 3 , and Mq. Taking m^ — ITeV, Mo = 
10 10 GeV, we plot the theoretical prediction of BR (T — y /i + r~) (solid line) 
and BR(r ^7) (dash line) versus 8f in FigfTJa). Taking M = 10 10 GeV, 
8f = 4 X 10 -3 , 8f = 3 X 10~ 4 , we plot the theoretical prediction of BR(T -» 
/U + r~) (solid line) and BR(r — y yrf) (dash line) versus in FigJT^b). The 
gray shadow is the excluded region for BR(r — y fij) from the experiment. 
Similarly, the mass insertion affects the theoretical evaluation of BR (T — y 
p + r~) strongly, and BR (T — y p + r~) depends on M mildly. 

The most stringent prediction of BR (T — y /i + r~) in theory is given in [lo[ 
with BR (T — y p + r~) < 2.9 x 10 -6 , and that is very close to the experiment 
limit. In [8], the author calculates the BR (T — y /U + r~) < 7.2 x 10 -5 , and it 
shows the updated data for BR(r — y 3/x) from experiment is necessary. Under 
our assumption, m^ = ITeV, M = 10 10 GeV, 5f = 4 x 10~ 3 , 5f = 3 x 10" 4 , 
we get BR (T —y p + r~) < 5.3 x 10~ 7 . That is also promising to detect 
directly in experiment in near future. 

We can evaluate the branching ratios of LFV decays p(u, J/^, T) —y e + p~~ 
using above method. To shorten the length of text, we just present the upper 
bounds on those branching ratios under the same assumptions as <fi — y e + pT . 
After considering the constraints from fj, — e conversion, /i —y ej, \i — y 3e etc, 
we give a summary of upper bounds of experiment data and corresponding 
theoretical predictions in TabJTJ 

4. Conclusions 

Considering the constraints from \x — e conversion, fi — y ej, r —y ^7, 
/U — y 3e etc, we analyze the LFV decays of — y e + /i~, J/^/(T) —y p + r~ in 
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the framework of MSSM with type I seesaw extended. 

In the MSSM with type I seesaw extended, the theoretical evaluation on 
BR (<fi e + /i _ ) is affected by the mass insertion Jf 3 ^ 3 . After considering 
the constraints from jjl — e conversion, \i — > and /i — > 3e, the prediction 
on BR (<f) — y e + fi~) can reach 5.0 x 10~ 20 , which is far below the present 
experimental upper bound. In a similar way, the mass insertion 5£ 3 affects 
the theoretical evaluations on BR (J/^ — » /i + T~) and BR (T — > 
sensitively. Considering the constraints from r — > /i7 and r — > 3/x, the 
predictions on BR (J/ V I / — > /i + T~) and BR (T — > h + t~) can reach 10~ 7 , 
which are very promising to be observed in near future experiment. 

In the future, the expected sensitivities for BR (ji — > ej) would be of 



order 10" 13 |35|. For BR (r ->- erf) and BR (r ->- /ij), it would be 10" 9 [36 
For CR (/i — e, Tz) in nuclei, it would be as low as 10~ 16 ~ 10~ 17 j37|. Thus, 
the n — e conversion experiments would represent the most promising tool 
to probe new physics. The study of LFV decays via vector mesons is also 
waiting for the new data from the experiment. 
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Appendix A. The simplified amplitude in seesaw extended MSSM 

In this appendix we present the simplified amplitudes in FigJTJ 

. o o P 6 4 



i,j=l k=l 



XU. 



+C 2 



(p,)\c ime (A kj A ik *P L + A^Af*P R ) 



m, 



,A k 3 j A\ k * - mil Af*A k /)P L 



+ (m e Af*A k ^ -m^Af^Pj 



R 



(m e Af*-m x oAf*)Al j P R 



+A k f(m e A^ - m xl Af*)P L ] }v e (p. 
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where 



Af = 



Aik 
A 2 



e(s w Z N + c w Z N 
c w s w 



2 rv3k rv'oi 



+ Y l Z N Z , 



Y?z?{zfy - 



y/2eZ*{Z$y 



k j _ e(s w Zj^ + CwZ^Z^ ! 3k 4j 
n 3 — /= r 1 1 Zj n Zj { 

= Y l 1 zf^(zl j y - ^^jH^lT. 



and 



with i — 1, 2. 



^ /^yr222 2 2 2 1 

Co = C [m ,m e ,m ,m L .,m L ,mJ 



2 _2 2 



^ r 2 2 2 2 2 2 

Ci = C i [m IM ,m e ,m <j) ,m Li ,m x o,m L 



where 



A, 



i&^f^m^sl - 3) 



(Pa + P4) • e(p) ^ 

j,j=i fc=i 



x4^ /i (p 4 ){c 1 m e (^Af *P L + A k jA?*P R ) 



+C 2 



m, 



+(m e A'l j Af* - m^A kj Af*)P R 



+C 



e^4 



+(m e A k 3 j Af * - m x o4Mf)? L ] } Ve (p 3 ) , 



2i nrlj* 



A l i = z?z; 



and 



n ^ r 2 2 2 2 2 2 l 

C = C [m ,m e ,m^m L .,m L .,m x o\ 

n n r 2 2 2 2 2 2 l 

Ci = Ci[m^ m e ,m^ m L .,m x o,m L .\ 
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A 



%e 7r 



7^(44 - 3) 



EE 



48 v / 2ct(m2-m2)iV cS t^ 1 ^ 
xuMfeip) [(m e m^A jk (Co + d + d) 
+2s w Af {m e m x o^(C + d) - m e m x o4d}W* 



+2s w A%*^ k {Ai i m x oC + (Agm^ - ^m x o)C 2 }m^ 
+Ai k {(4A$ - 2Ai i m x om x o)C - A$s w [l + 2C 12 mJ 



-2(Ci + Ci 2 + Cn)mJ - 2(d + C 12 + d^m*]} 
+2(p 4 • e(p)) [(A* A* {(d + di)m e P L + Ci 2 m M Pfl} 
+2A^f Sw m x odPii)4 fe * - 2s w (4 i 4 fc m x oC 1 P L 
+2^Af Sw {d 2 m M P L + (d + di)m e P R })Af *" 
+2(p 3 • e(p)) [2^ Sw (4 fe m x odPL + 2A{ k s w {(d 
+C 22 )m,P L + Ci 2 m c P fl }) A ik * - Ag (A?C 12 m e P L 



+{Af(d + C 22 )m^ + 2A| fc Sw m x od} J P«)^ fc * 



-^(p) [{AgA? [(d + d 2 + d 2 K - C 12 mJ + (d 
+di + d 2 )mj]) + 2 ({24* - A%}m x om x oC 
+AAi k {A^m x o(C + d) - ^m^ 2 }m^ w )4* 
+4m e s w (4'Af m x oC + A{ k (A^m x o - Agm^d 
+2A«Ai k (C + d + C 2 )m„s w )A[ k * - l] P L }^(p 3 ) 



^6 — N N 



A jk = Zl k (Z 



TV °w 



z# Cw ) + r/zf z 



A' 



^8 N + J Z L N 
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\ik r?5k iyli* 1 -1^2 ^2fc ?73i* 

^10 — + r « N ■ 



and 



C = C [m^, m e , m M , m Xi , m^o, m L J 
Ci,ij = Ci^ [mj, mj, ml, m\ k , m 2 , m 



with i — 1, 2. 



.A 



where 



ie 2 7r 2 /^ 



EE 



w ij=i fe=i 



x^(p 4 ){^(p) [4i*({2(C 2 + C 22 )m\ + 2(d + C u )m 
+2C 12 (m 2 + mj - mj) + 2C (m x ±m x ± + 2) - l}e 2 Ag 
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A[ k = Y^Z 2 T(Z 2k - IZf) 



and 



A{1 = Z!{(Z?-IZ?) 
A* = Y?Z*\Z?-IZ?) 



r , ^ r 2 2 2 2 2 2 1 

ri ^ r 2 2 2 2 2 2 l 



with i — 1,2. 

Ae - 



m 2 e 2 f ^m^sj - 3) 
48c* (m| - ml)N c si 



EE 

i,j=l k=l 
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12 



mi 
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+4(p 3 • e(p)) eA***Ay 6 ( eA{*C 12 m e P L + {e^(C 2 + C 22 )m M 
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-s w A{ k 4 C 2 m x ±}P R ) + s w A$A¥ B \8 w A{i{(C 2 + C 22 )m,P L 
+C 12 m e P R } - eA{ k 3 C 2 m x ±P L ^ }v e (p 3 ) 



where 
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n /-< r 2 2 2 2 2 2 l 
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ri ^ r 2 2 2 2 2 2 l 
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C = C7o[m^,m e ,m j[1 ,m I/ .,m 1/j .,m x ±J 

/"< ri r 2 2 2 2 2 2 l 



1,2 



the simplified amplitudes in FigJ2J 



•An 



48JV,s; 



E EE 

i,m=l j=l 1=1 



*u,{pM2UmAA™ieD A^A%{em e m^At 



25 



where 



and 



+s w m e m xt Ay 2 *) - A^ls^A^A%{em xt m^l[ 
+ s ^ m xt m x^ A i2) ~ KislAtA^A%(2C, 
+2D 1 ml i )]t{p)PR + 2 i>4 - A™le 2 AtA™fA% 
x(2C + 2/W ) + AZ j D A^*Af 8 {s w (em e m x ± 
xAY* + m e m^s w A^) - e(m ±m M s w A l / 2 * + em ± 



Xm 

xm^A\{*)} {( P )PL + D fjA™ l *A% eA^{ 



em e 

xm x ±A\\* + m e m^s w A\^) - A™1 s w (em x ±m x ±A\ J * 
+m x ±m,s w Ag*)] {4{p)i> - Mp))Pr + A)/J<*4< 7 
x s w A™l (em e m^A l {* + m e m x ±s w A l / 2 *) - eA 7 ^ 
x (em x ±m^A\\* + m x ±m x ±s w A\ J *) tf(p)$ - #(p)) 
•^l}^(p 3 ) 



/i 17 

^18 



6 rvJU rvl(i,m) % .., y 

6 W 



J v (J+3)l ^2{i,m) 



K 



.12 



y2 gjl* r?2(i,m)* 



K 



.12 



n r 2 2 2 2 2 2 i 

Co = C^m^m^m^m^.m^.m^^ 



n r 2 2 2 2 2 2 2 2 2 2l 

D [m e , m s , m s ,m^, m D) m^ m Vj) m x ±,m Ui) m x ±\ 

n nr2 222222 2 2 2i 
L>1 = DoFd, ™ S . m e. m S , %± . m ^±l 



where m 2 D = m 2 ^ + m 2 s — 0.5m 2 



ITX 



48 JV, 



EE 



x^(j9 4 )|2/ m^Ar(<* J Do<m x o(^m e 
-Ar'm x o ) - <Mf (2C Ar + 2/W'mL)) 



26 
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